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We present a numerical study of the process of the kink-antikink collisions in the coupled one- 
dimensional two-component (f) 4 model. Our results reveal two different soliton solutions which 
represent double kink configuration and kink-non-topological soliton (lump) bound state. Collision 
of these solitons leads to very reach resonance structure which is related to reversible energy exchange 
between the kinks, non-topological solitons and the internal vibrational modes. Various channels 
of the collisions are discussed, it is shown there is a new type of self-similar fractal structure which 
appears in the collisions of the relativistic kinks, there the width of the resonance windows increases 
with the increase of the impact velocity. An analytical approximation scheme is discussed in the limit 
of the perturbative coupling between the sectors. Considering the spectrum of linear fluctuations 
around the solitons we found that the double kink configuration is unstable if the coupling constant 
between the sectors is negative. 

I. INTRODUCTION 

Since the early 1960s, the soliton solutions in the non-linear field theories have been intensively studied in various 
frameworks. It is evident these spatially localized non-perturbative configurations play important role in a wide 
variety of physical systems. The study of the interaction between the solitons and their dynamical properties has 
attracted a lot of attention in many different contexts. 

Simplest example of the topological solitons in one dimension is the class of the kink (K) solutions which appears 
in the model with a potential with two or more degenerated minima. Double well potential corresponds to the nonin- 
tegrable cf) A model. Further extentions of this model are possible, for example one can consider domain walls obtained 
by embedding of the 1+1-dimensional <j> -kink into higher dimensions. This model has a number of applications in 
condensed matter physics [T] , field theory [3] and cosmology [J] . Probably one of the most fascinating developments 
in this direction is the idea that one may regard our universe as a domain wall [3 [6] and Big Bang is associated with 
the collision of the domain walls [7]. 

Further, in various branches of physics more complex systems involving several coupled scalar fields arise (see, e.q. 
[8HTU]). Properties of the solitons in these extended models are rather different from being a trivial extension of 
the single-component one, in particular the model may support existence of non-topological solitons [8] or it can be 
related with application of the supersymmetry in nonrelativistic quantum mechanics However, relatively little 
is known about the dynamics of the multicomponent configurations in the extended non-integrable models. 

Dynamical properties of the usual 4 kinks, the processes of their scattering, radiation and annihilation have 
already been discussed in a number of papers, see e.g. [12H19j . Surprisely, it was discovered numerically that there 
is dynamical fractal structure in the process of the kink-antikink (KK) collision P31 HH HE1 HH] • Depending on the 
impact velocity, the collision may produce various results, for some range of values of the impact velocity the KK 
pair form an oscillating n-bouncing window, after collision the pair reflects back to some finite separation and returns 
to collide (n — 1) times more before escaping to infinity. Also the process of production of kink-antikink pairs in the 
collision of particle- like states is chaotic [Mj. Another interesting phenomenon occurs when the kink configuration is 
under the influence of an incident wave |21) , it was found that the <f> kink starts to accelerate towards the incoming 
wave although the effect depends on the structure of the potential. 

The explanation of the appearance of the resonant scattering windows in the KK collision is related to the energy 
exchange between the translational collective mode of the solitons and an internal vibrational mode of the kink 
[HI H3 HFJ [TS]. Another mechanism works in the cf> 6 model where the resonance windows appear due to energy 
transfer between collective bound states trapped by the KK pair [22) . Evidently, in the extended multicomponent 
model coupling to the second component may strongly affect the process of the resonant energy transfer. 

A natural extension of the simple model </> 4 model in 1+1 dimension is related to consideration of systems with 
two or more component scalar fields. Several models of that type have been investigated in recent years, some of 
them are related to consideration of the complex scalar field, a deformation of the linear 0(2)-sigma model known 
as MSTB-modcl 23 25 or with consideration of a two-component non- linear interaction model [51 151 12"6l - f2"8] . In 
particular, this model supports existence of solutions with one of the component having the kink structure and the 
second component being a non-topological soliton. Another situation occurs when the system of coupled scalar fields 
belongs to the bosonic sector of a supersymmetric system (Wess-Zumino model with polynomic superpotential and 
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two Majorana spinor fields) or with application of supersymmetric methods [TT] and related reduction of the 
second-order field equations to the system of corresponding first-order Bogomol'nyi equations [30 | 131 ] . 

The aim of the present paper is to analyse the process of the KK resonant bouncing scattering in the two-component 
coupled 4 model. We assume that both sectors are two replica of the usual 4 model, the fields are coupled through 
the minimal quadratic term and the kink-antikink collisions are taking place in the first sector. Evidently the fine 
mechanism of the reversible energy transfer in the bouncing of the kink will be modified in such a model because both 
the vacuum structure and the spectrum of linear oscillations around the solitons will be different from the standard 
one-component theory. 

This present work is organized in the following way. In Section II we overview the spectrum of soliton solutions of 
the coupled two-component 4 model and discuss the analytical perturbative approximation to these configurations in 
both sectors. In section III we present our numerical results of the various kink-antikink (KK) collisions and discuss 
the resonance structures we observed. In section IV we describe the spectral structure of linear pertubations on the 
background of the soliton solutions of this model and discuss the possible mechanisms of the resonance energy transfer 
between the internal modes and the kinks. We present our conclusions in section V. 

II. TWO-COMPONENT COUPLED </> 4 MODEL 

We consider the system of two coupled copies of the 4 fields 

L = \d^id^! + \d^ 2 d^ 2 -\{4>\- l) 2 - \ (4> 2 2 - a 2 f + k&<& . (1) 

where a is associated with the mass of the field 2 and n is the coupling parameter |38j . Thus, the potential of the 
model is 

U(<h,<M = ~ (0? ~ if + ~ (0l ~ a 2 ) 2 - (2) 

where k is the coupling constant. Hereafter, for the sake of simplicity we set a = 1. 
Let us consider the vacuum manifold of the model which is defined by the conditions 

01,* = <h,t = <f>i, x = 4>%x = 0, £/(0i,0 2 ) = min, (3) 

that is 

20i(0?-l-/^) = O, (4a) 

202(0i - l-K^i) =0- (4b) 

Evidently, the structure of the vacuum depends on the value of the coupling constant n. In the case when — 1 < k < 1 
there are nine stationary points of the potential ([2]): 

• one local maximum 0i = 02 = (U max = 1), 

• four saddle points 0i = 0, 2 = ±1 or <j> x = ±1, 2 = {U sad die = |), 

• four minima 0i — ±l/\/l — k, 02 = il/yl — n (U m i n = — jz^t)- 

If k > 1 the system becomes unbounded from below, that is it is unstable, so we shall not consider this case. In the 
limiting case k = — 1 the minima of the potential |2]) and its saddle points are degenerated while in the regime n < — 1 
the minima appear at the four points 4>\ = O,0 2 = ±1 or <f>\ = ±1,02 = 0, which previously, when —1 < k < 1, 
were associated with the saddle points. Correspondingly, the former minima become the saddle points in this case, 
see Fig[T] 

Thus, there are kink soliton solutions in this model which interpolate between one of the vacua at x = — oo and 
another one at x — +oo. To find these configurations explicitly we have to solve the Euler-Lagrange equations of 
motion, which are the following coupled, nonlinear, partial differential equations: 



d tt (pi ~ 3 ra 01 + 20! (0? - 1 - K 02) = 0, 
dttfa ~ 3 X ,02 + 202(0^ - 1 - K0f) = 0. 



(5) 
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Figure 1: (Color online) Contour plot of the potential (J2J for k — 0.5 < 1 (left panel) and k = —1.1 < —1 (right panel). 



Let us consider |, 

(±i/vT=«,±i/vr 

topological. However the set of discrete symmetries of the Lagrangian ([I]): 



< 1. Then a particular static kink solutions interpolate between two of four vacua: 
c). So there are 2( 4 ) = 12 types of the (anti)kinks in the model, both topological and non- 
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actually reduces the number of significantly different solitons to two. 

Indeed, the solitons can be classified according to their topological charge 



dx 



dx 



1,2. 
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For each component it can have three values {—1,0,1}, kink soliton has topological charge +1 while antikink has 
topological charge —1. Thus, the configuration can be labeled by the topological charges of the components as (n\, 712) 
where m £ {—1, 0, 1}. 

Let us consider the most symmetric case 4>2 — ±01 = 4>- Then the system of coupled equations ^ becomes reduced 
to equation 



d U (f>- d xx (f> + 2(/)((l - k) 



1) 



0. 



(8) 



Evidently, rescaling of the field cf> — >■ (f> = <j> y/l — k yields the familiar equation of motion of the usual one-component 
4 theory. So the model |TJ) in this case is reduced to the two separable copies of the 4 model. Static solutions 
{dt4> = 0) of the system ^ can therefore be written as 
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(9) 



k as x increases from 
k and — 1 / \/l 



00 to 00 
k as x varies from 



The static kink solution interpolates between the vacua — k and l/\/T 

Similarly, the antikink solution interpolates between the stable vacua 
to 00. The corresponding 2-component configuration possess non-zero topological charge in both sectors, so there 
are four kinks which correspond to the configurations with topological charges (1, 1), (—1, 1), (1, —1) and (—1, — 1), 
respectively. We shall refer to this solution as double kink. 

The perturbative sector of the model consists of small linear perturbations around one of the solutions, thus by 
analogy with the usual one-component 4 model, one can expect the process of collision between the solitons is related 
to reversible exchange of energy between the translational and vibrational modes of the individual kinks [Til [Pol [TBI [T§] . 
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Figure 2: Lump kink static configurations (1,0) is plotted for k = 0.5. 



However in our case there is another possibility related to energy exchange between the sectors, so the dynamics of 
the kinks and the spectral structure of the perturbations can be rather interesting. 

Besides the double kink configuration, there is another solution to the model ([I]) which corresponds to the topological 
soliton in one of the sectors and non-topological soliton (lump) in another sector. We shall refer to this solution as 
lump kinks. In such two-component system the field of the topological soliton in one of the sectors interpolates 
between the vacua =f1/>/1 — k and ±l/\/l — n as before, while the non-topological soliton belongs to the sector with 
topological charge 0, each of these lump solitons can be settled in one of two vacuum states, so there are 8 such 
configurations labeled by the topological charges (1,0), (—1,0), (0, 1) and (0, —1). For example, the static lump kink 
solution having first component interpolating between the vacua — n and 1/yl — k and the second component 

being the non-topological soliton about one of the vacua ±1 /a/1 — K is presented in Fig. [2] 

Other static lump (anti)kinks can be obtained from this solution by using the discrete symmetries of the model 
([6]). Qualitatively, one can understand the physical reason of the lump kink solitons existence, if we consider the kink 
configuration in one of the sectors, in the configuration space (4>i, fa) it tends to evolve from one vacuum to another 
through the saddle point (fa(0) = 0,^(0) = 1). Thus, non-topological soliton in the second sector appears. However 
the resulting configuration does not reach this point exactly because of the field derivatives contributions in the total 
energy. Note that although the general structure of the solution is similar to the configuration considered in [8| , there 
is no analytical solution of the model ([I]), the kink component there is rather different from the solution of the usual 
fa model. In the case under consideration the topological solitons can be treated as "stretched out" fa kinks which 
are coupled to the lumps in the second sector. These configurations can be constructed numerically as solutions of 
the of static equations ^ when we impose the boundary conditions 

fa(0) = 0, 0i (00) = 1/vT^, 4(0) = 0, fa(oo) = i/vT=«. (10) 

Since the coupling constant |«| < 1, we can construct an analytical approximation of the lump kink solution 
considering the coupling between the components perturbatively. The initial decoupled static configuration consist of 
a kink in the first sector of the model ([lj and the trivial vacuum solution in the second sector: 

=tanh(a;); 4° } = 1- ( n ) 

Corrections to this configuration can be obtained from the scries expansion of the fields fa , fa in powers of perturbation 
constant [21] [33] 



fa = + K 0« + + . . . . ^ = (O) + K0 (l) + k 2 (2) + 



(12) 
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Figure 3: The mass of the two kink types, the lump kink and the double kink, as function of coupling constant k. 



The first-order time-independent corrections to the initial solution (111 then can be obtained from the equations 

+ (6tanh 2 x - 2) <^ 1} - 2 tanh x = 0, 



U^) + A^p ~ 2 tanh2 x = °- 

V / XX 



(13) 



Evidently, as x — > oo the first order correction (jfp — ¥ | supposing that on the spacial asymptotic Ui j 0. This 

V / XX 

is in agreement with the vacuum shift due to perturbative coupling in both sectors — n ~ 1 + |ft. 

Further, we can try to construct an analytical approximation to the first order correction using an expansion in 
hyperbolic functions 



J2 A n tanh 2n+1 . 



E 
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cosh 2 ™ x 



(14) 



Substitution of these series into the first order equations ( 13 ) yields 

4 



!>$ 1 ^ = — tanh a; tanh 3 x - 



16 

15 
1 



75 



tanh 5 



x tanh 7 x 

75 



1 



1 



45 



(15) 



2 2 6 cosh 2 a; 28 cosh 4 a; 728 cosh 6 x 11284 cosh 8 a; 
Similar approximation can be applied to the higher order corrections. 
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Note that the perturbative expansion reveals instability of the double kink configuration in the case k < 0, the 
system may decay into a pair of lump kinks. Indeed, the Figure [3] shows the masses of kinks as a function of coupling 
constant k. Evidently, for negative values of the coupling k the mass of a double kink is less than the mass of two lump 
kinks. That means that the double kink is less energetically favorable state than two separated lump kinks. Kinks 
in one-component 4 model are stable, therefore the perturbation which can destroy the double kink configuration 
must break the symmetry <f>\ — ±<fo- For positive k the double kink is energetically stable. Moreover it is possible 
that the double kink could be a final state in collisions of two lump kinks with kinks in different sectors, the outcome 
evidently depends on the parameters of the model, the impact velocity and the coupling k. Below we discuss various 
solitons collisions scenarios and classify them. 
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III. KINK-ANTIKINK COLLISIONS IN TWO-COMPONENT <j> 4 MODEL 

A phenomenon known as the two-bounce resonance soliton scattering in the usual one-dimensional </> model was 
observed by many researchers in numerical simulations [H] [HI HH [19] . Unfortunately there is no complete analytic 
description for such process, the model is non-integrable and the dynamics of the solitons is highly complicated. 
Numerical simulations performed over last 20 years [21 HU HH1 HH] discovered an interesting result: the process of 
collision between the solitons is of fractal nature, depending on the impact velocity, it may produce various results 
and for some range of values of the impact velocity the KK pair form a set of bouncing windows. 

Although semiclassical consideration [TSl [TH] [H] could capture the most important features of the process, the 
analysis entirely rely on extensive numerical computations. It is known that for initial velocities of a kink and 
an antikink above a critical value v cr — 0.2598 two incident solitons always escape to infinity after collision, with 
the emission of some radiation, thus the process is similar to quasielastic scattering of two particle-like objects 
[21 [T51 [THl HH] ■ For the impact velocities below v cr , the kink-antikink pair generically becomes trapped, it annihilates 
into radiation, however numerical analysis reveals the fractal structure with a sequence of narrow bouncing windows 
within which the solitons are again able to escape to infinity. The first of two-bounce window opens at the impact 
velocity w in = 0.189 then the resonance windows become more narrow as the impact velocity increases. 

If we look on the edge of a two-bounce window, we can observe a narrow sequence of the 3-bounce windows which 
are characterized by the same relations between the windows widths and the impact velocities (self-similar structure) 
[18j . These windows are called 3-bounce ones. Furthermore, on the edge of every 3-bounce window we can find the 
4-bounce windows structure, etc. However, the fractal dimension of this dynamical structure of the collision of the 
solitons is restricted, so we refer to this to as a quasi-fractal structure. 

The accepted explanation of the appearance of these windows is that they are related to reversible exchange of 
energy between the translational and vibrational modes of the individual kinks. At the initial impact, some kinetic 
energy is transferred into internal shape modes of the kink and antikink. They then separate and propagate almost 
independently, but for initial velocities less than v cr they no longer have enough translational energy to escape their 
mutual attraction, and so they return and collide one more time [HI [18] . At this point some of the energy stored 
in the shape modes can be returned to the translational modes, tipping the energy balance back again and allowing 
the kink and antikink to escape to infinity, provided that there is an appropriate resonance between the time interval 
between the two collisions, and the period of the internal modes. More generally, sufficient energy might be returned 
to the translational modes after three or more kink-antikink collisions, leading to an intricate nested structure of 
resonance windows. Thus, the resonance condition is that the time T between two consequent collisions should satisfy 
the equation [TS] [TH] 

lo s T = 8 + 2i:n (16) 



where 5 is a phase, and n is an integer [39]. It can be shown that ^ = 5.13 and S = 3.3(« n). The relation (16) can 
be modulated [57] , 

This mechanism allows us to construct an effective model of the collision process [15]. The idea is to separate 
collective coordinates of the kink-antikink system and make use of the Yukawa interaction between the well separated 
pair mediated by the meson state of mass m. Surprisingly enough, this model works rather well. 

Our goal is to analyse the process of collision of the soliton configurations in two-component model (JlJ numerically. 
The initial data used in our simulations represent widely separated kinks or lump kinks propagating from both sides 
towards a collision point, the parameters of the collision process are the impact velocity v and the coupling constant k. 
To find a numerical solution of the PDE describing the evolution of the system, we used the pseudo-spectral method 
on a discrete grid containing 2048 nodes with periodic boundary conditions. For the time stepping function we used 
symplectic (or geometric) integrator of 8th order to ensure that the energy is conserved. The time and the spacial 
steps are 5x = 0.25 and St = 0.025, respectively, so the numerical errors are of order of (St) s . 

Evidently, if the sectors of the two-component model are decoupled, the pattern of collision reproduces the 
chaotic structure we described above. 



A. Lump kink collisions - same sectors 



Let us consider collision of two lump kinks (—1,0) + (1,0) at k = 0.5. The numerical simulations reveal rather 
interesting pattern, for initial velocities of a lump kink and an lump antikink above a critical value v cr — 0.453, an 
intermediate state which is created in both sectors at the collision center, decays into lumps in the first sector and 
into kinks in the second one. Then two solitons escape to infinity, with the emission of some radiation, so the process 
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Figure 4: Collision of two lump kinks (—1, 0) + (1, 0) — > (0, — 1) + (0, 1) at « = 0.5 and initial velocity Vi n = 0.5. 



is similar to quasielastic scattering of two 4 kinks up to the flip between sectors: 

(-1,0) + (1,0) ^(0,-1) + (0,1). 

This process is displayed at Fig. [4j For the impact velocities below v cr we the observe the resonance energy exchange 
between the sectors with a sequences of bouncing windows within which the solitons are again able to escape to 
infinity, as seen in Fig. [5j The first of two-bounce windows open at the impact velocity Vi n = 0.423, the consequent 
structure of the resonance windows is not so regular as in the case of the usual 4 model, also some false windows 
appear. Evidently, this can be related to more complicated structure of the spectrum of excitations, alongside with the 
internal mode of the kinks there are degrees of freedom associated with excitations of the lumps, also some collective 
degrees of freedom in the system kink-lump may appear. 

In all windows we observed the same permutation of the kink and the lump sectors in the final state. For smaller 
impact velocities and between the windows, the intermediate state which is created in both sectors at the collision 
center, decays into two oscillons (one in each sector) with some amount of radiation emitted. The final states of the 
process then are: 

! (1,0) + (-1,0) bounce (I) 
(0,1) + (0,-1) flip (II) (17) 
(0, 0) annihilation into oscillon (III) . 

However for ultra-relativistic impact velocities (for example at v = 0.950, n = 0.386, see Fig. [8]), we also found another 
final state 

(1,0) + (-1,0) ^(-1,-1) + (1,1) (IV). (18) 

In this process relativistic collision of the lumps leads to the production of the kinks in the second sector. They bind 
to the kinks from the first sector and then the pairs travel as double kinks. 

In Fig. [7] we present the final vacuum state in the center of collision after a long time 4n{x = 0, t = 150) as a 
function of the parameters of the process of collision (1,0) + (—1,0), the value of the coupling constant k and the 
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Figure 5: (Color online) Two lump kinks (—1,0) + (1,0) collision at k = 0.5. The plots represent the field values measured 
at the collision center as function of impact velocity and time. Upper and bottom plots represent the <f>\ and 02 components, 
respectively. 

initial velocity Vi. Before the collision the initial value of the field at the center was taken as 0i(O, 0) = -^== - (the 

only possiblity for such configuration) and 02(0,0) = — ^zr^ ( one °f two possibilities). We will denote the state of 

that vacuum in square brackets as $ = \A — K [<^i; 4>2\ = [1, — 1]- The final state of the vacuum in the collision center 
then can be defined as $ = [— 1, — 1] + 2Q, where Q is the topological charge of the kink moving towards asymptotic 
x — > — co. 

Depending on the parameters of the collision process, we observed various scenario. If the final vacuum state 
remains the same after collision, the kinks are bounced back (cf. Fig. |8j panel (I) , where we plotted positions of the 
topological zeros for that collision). On the upper left plot in Fig. [7] the bouncing corresponds to the regions colored 
in blue. The largest of these regions has a shape of a triangle at k = 0, note that the bouncing is also observed for 
some set of higher impact velocities and relatively strong values of the coupling k, both positive and negative. 

Black regions in Fig. [7] correspond to the final state [—1,1], in other words the kink propagating on the left 
towards the spaciai asymptotic x —> — oo after the collision now belongs to the sector with topological charge Q = 
|([— 1, 1] — [— 1. —1]) = (0, 1). Thus, the collision leads to the topological fliping between the sectors (cf. Fig. ^1 panel 
(II)). 

The regions in Fig. [7] which are colored in red, correspond to final state of the configuration in the region between 
the collision center and the asymptotic x — > — oo with total topological charge equal zero. This is the process of KK 
annihilation which dominates at relatively small initial velocities (cf. Fig. [8j panel (III)). 

Finally, we observed that for the relativistic impact velocities and relatively large positive values of the coupling 
constant k the process of the production of the double kink is allowed as illustrated in Fig. [8j panel (IV). The 
corresponding regions belong to the sector with topological charge Q = (1, 1), in Fig. [7]they are colored in yellow. 

Note that the borders between those regions are very complicated, presumably they are demonstrating some kind 
of fractal structure. 



B. Lump kink collisions - different sectors 

Since the kink is complemented by a non-topological partner in the opposite sector, there is another process of 
collision between the kinks from different sectors, mediated by the lumps. Numerical simulations reveal that at 
relatively small velocities v < 0.565 the energy of collision between the kink and the lump is transformed into 
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Figure 6: (Color online) Collision of two lump kinks in different sectors mediated by the lumps at k = 0.5. 




Figure 7: (Color online) Final state of the two lump kinks collision (1,0) + ( — 1,0). 
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Figure 8: (Color online) Example of topological zeros motion for two lump kinks collision (1,0) + (—1,0) in the same sectors. 
(I) bounce (1,0) + (-1, 0), (II) flip (0, -1) + (0, 1), (III) annihilation (0,0), (IV) double kinks creation (-1, -1) + (1, 1). 



excitations of the kinks which are erratically oscillating around the collision center it the both sectors, see Fig. [6] 
Thus the resulting configuration is actually an excited double kink. It will be shown in the next section that the 
double kink for k > has at least two oscillational modes. One of them corresponds to the oscillational mode of the 
<p 4 kink and has the same frequency \/3. There is also a mode which is responsible for translational oscillations of the 
kinks in both sectors around the center of mass. Both of these frequencies can be seen in the power spectrum of the 
field measured in the center of mass. Moreover due to nonlinearities combinations of those frequencies are also visible 
(cf. Fig.§). 

For impact velocity v — 0.565 we observed the bouncing of the kinks, a narrow resonance window opens, a sequence 
of the bouncing windows of decreasing width follows it up to the v cr = 0.584. This resonance structure is much 
more regular than in the case of the collision of two lump kinks in the same sector, there is no false windows and 
the dependence of the windows on the impact velocity closely reminds the well known pattern observed in the usual 
one-component 4> A theory |15l I18j . Evidently we can consider it as undirect evidence of the simple mechanism of the 
resonance energy exchange between the kink and the lump which involves only the internal mode of the kink, there 
is no energy exchange between sectors. 

Above this threshold velocity the lump kinks are passing through the lumps associated with the kinks from the 
second sector, thus this is a regime of quasi-elastic collision, see Fig. [T0| 

Variations of the parameters of the collision process, the coupling constant n and the impact velocity v reveal several 
possible final states: 

!(1, — 1) double kink creation (I) 

(0,-1) + (1,0) passage (II) (19) 
(1,0) + (0,-1) bounce (III) 

Note that the double kink production is allowed only for k > 0, we already mentioned that for k < the double kinks 
are unstable with respect to decay into the lump kinks. 

All the processes decribed above are illustrated in Fig. [TT| which presents the final state of the collision (1,0) + (0, 1) 
associated to the kink on the left hand side from the collision center propagating towards the spacial asymptotic 
x — > — oo as $ = [—1, —1] + 2Q. Similar to the Fig. [7J the black regions $ = [—1, 1] correspond to the quasi-elastic 
scattering of the kinks (passage), in this case the topological charge of the scattered kink in the final state is Q = (1, 0) 



(cf. Fig. 11 (II)). Then the regions in Fig. 11 which are colored in blue, correspond to the bouncing of the kinks, then 



in the final state we have Q = (1,0) (cf. Fig. 11, regions (III) and (IV)). Note that there is another possibility, as a 
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Figure 9: (Color online) Power spectrum of the field in the first sector in the collision center for (1,0) + (0,-1) — > (1,-1) 
and k = 0.5 at Vi n = 0.2. As the result of the collision an excited double kink is created. One can identify the frequency of 



the usual 4 internal mode uj^i = \/3, another frequency 



1.129 (see the next section for details) is 



lowered because of nonlinearities to u)2 = 1.077, other combination of those frequencies. 



result of the collision two kinks can merge forming an (excited) double kink. These processes are depicted as whitish, 
yellow and light blue colored areas on the right side of the plot (Fig. 11 (I)). 



C. Double kink - lump kink collisions 



Another interesting solitonic collision process in the two-component model ([!]) is the double kink - lump kink 
collision. Let us take the initial state as, for example (—1,-1) + (1,0). Since the total topological charge in the 
first sector is we could expect that the outcome of the collision, at least at small impact velocities, will be just 
annihilation of the lump kink and one of the components of the double kink. However, a single component of the 
double kink is not a solution of the model Q, so it has to deform into the lump kink then. Indeed, at relatively 
small impact velocities v < 0.76 we observed such a process (— 1, —1) + (1, 0) — > (0, —1) + (0, 0), besides radiation and 
oscillon remains in the first sector as the result of annihilation of the kink-antikink components. 

It is known that the resonance excitation of the oscillon can produce KK pairs [20], so at higher energies other 
scenarios are possible. Indeed at K = 0.5, within the impact velocity range 0.76 < v < 0.82 we observed another 
resonance structure between 3 channels 



(-1,-1) + (1,0) 




(0,0) 

(0,1) + (0,-1) 
(1,0) + (-1,0) 



(20) 



The results of numerical simulations are presented in Fig. [13j where we plotted final velocity of the solitons vs the 
impact velocity. The bottom plot represents zoom in on the smaller subrange of the impact velocities in between, 
evidently this plot is qualitatively similar to the upper plot. Thus a new type of self-similar fractal structure appears. 
Note the width of the windows increases with the increasing of the impact velocity, so this structure is inverted with 
respect to the usual fractal structure observed in the (j) 4 model [H] [18] . The complicated fractal structure probably 
is possible due to the fact that the double kink has at least two bound modes which can couple to each other and to 
the oscillational mode of the lump kink. Since at the collision process the double kink behaves as a bound state of 
the lump kinks, this type of collision can be interpreted as a three body problem. 
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Figure 10: (Color online) Collision of two lump kinks in different sectors mediated by the lumps at k = 0.5. The plots represent 
the field values measured at the collision center as function of impact velocity and time. Upper and bottom plots represent the 
4>i and <f>2 components, respectively. 




Figure 11: (Color online) Final state of two lump kinks collision (1,0) + (0, 1) in different sectors. 



Note that for k < the double kink is unstable, however the perturbation due to the lump kink is relatively weak, 
so the double kink is moving as a single object and it decays into the pair of lump kinks just before or already during 
the collision. 

Finally, as the impact velocity increases above the critical value v cr = 0.82, quasi-elastic scattering with flipping of 
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Figure 12: (Color online) Example of topological zeros motion of two lump kinks collision (1,0) + (0, 1) in different sectors. (I) 
double kink creation (1, 1), (II) passage (0, 1) + (1, 0), (III) bounce (1, 0) + (0, 1), (IV) multibounce (1, 0) + (0, 1). 




Figure 13: (Color online) Double kink - lump kink (—1, —1) + (1, 0) collision at k = 0.5. The plots represent the escape velocity 
of different solitons produced in both sectors after collision as a function of impact velocity. Bottom plot represent the zoomed 
in region of the upper plot. 
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Figure 14: (Color online) Final state of the double kink with a lump kink collision (1, 1) + ( — 1, 0). 



the sectors is observed: 



(-1,-1) + (1,0) ^(+1,-1) + (-1,0) 



(21) 



thus, the total topological charge in each sector remains the same but the kink component of the lump kink is getting 
captured into the double-kink state while the second component of the former configuration becomes the lump kink. 

lump kink collison (1, 



In Fig. 
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: we have gathered our results for the various channels of the double kink - lump kink collison (1, l)+(— 1, 0). 
Here we present the values of the fields at the center of collision, evidently the spectrum of the final states is much 
richer than in the case of the lump kink collisions we discussed above. The colors which represent the final state in 
Fig. [l4j are similar to the palette we used in the Figs. |7|lT] above. Note that since the double kink and the lump kink 
has different masses, the center of collision is not a center of mass of the system, so the dynamics becomes rather 
involved. For the sake of completeness in Fig. 15 we also presented various examples of the evolution of the topological 
zeros, associated with the position of the solitons in the process of the double kink-lump kink collision. 



IV. SPECTRAL STRUCTURE OF LINEAR PERTURBATIONS 



There is certain similarity between chaotic structures observed in the collision of the kinks in the two-component 
model |l]) and KK resonance bouncing in the usual 4 model [TS1 [TS] . This analogy can be evidently seen in the 
collision of the lump kinks, however the chaotic structure of the collision between the double kink and the lump kink 
is much more complicated. 

To investigate this analogy in more details we have to make use of the linear expansion around the solitons of the 
model under consideration assuming perturbations of both components (4>i +?7ie 1 "', </>2 + f?2e iti; *), here the unperturbed 
configuration (</>i, fa) corresponds to one of the solutions of the model ([I]), the double kink or the lump kink. 

Linearization of the system ^ yields the matrix eigenvalues equation 



-Ak4> 



192 





m 


= oj 2 






r i 2 . 




P 2 . 



where the operators D\ , D\ are 



Dl 



dx 2 
d 2 
dx 2 



b 2 — 2 — 2k</>2, 



Hi 



2 - 2k, 



(22) 

(23a) 
(23b) 
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Figure 15: (Color online) Example of topological zeros motion for collision of a double kink with a lump kink (1, 1) + (—1,0) 
with final states (I) unstable double kink decays to pair of lump kinks, one is ejected and second annihilates with the other kink 
(1, 1) + ( — 1, 0) — > (0, 1) + (1, 0) + (— 1, 0) — > (0, 1) + (0, 0), (II) annihilation in the first sector and passage in the second sector, 

(III) annihilation in the first sector, the created oscillon is coupled to the remaining kink in the second sector (0,0) + (0, 1), 

(IV) bounce (1, 1) + (—1, 0), (V) annihilation in the first sector and bounce in the second sector (0, 1) + (0, 0), (VI) bounce in 
one sector and passage in the second sector (1, 0) + (— 1, 1). 



Note that in this section we do not consider the coupling constant n as a perturbation parameter. 



A. Spectral structure of double kink 

First, let us consider the spectrum of linear oscillations around the double kink solution. In this case <p2 = = 



± 



tanh x 



and it is possible to diagonalise the system ( 22 ) introducing new functions 



£2 = 2 fa - 



(24) 



Here the variable £1 describes the motion of the center of mass of the system whereas the variable £2 corresponds to 
the oscillations of the components about it. Then the system of equations (22) is decoupled and we get 



dx 2 
d 2 6 + 2k 
dx 2 1 — K 



6 tanh 2 x - 2 £1 = ^ 2 £,i 



tanh 2 x - 2 f 2 = w 2 £> 



(25a) 
(25b) 
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Evidently, the first of these equations coincides with the usual modified Peashle- Teller equation that describes the 
linear excitations around the (j) 4 kink, it is well known the spectrum of fluctuations in this case has one translational 
zero mode cj = 0, the internal mode of the kink lu 2 = 3, which corresponds to the oscillations of the width of the kink, 
and the continuum modes ui 2 = k 2 + 4: 



(o) _ 



cosh 



(i) _ 



sinh : 



si 



cosh x 

(3tanh 2 x — "iiktanhx — 1 — k 2 ) 



(26) 



vi 0) 



The second linearly independent solution is the complex conjugation of the above. 

In the case under consideration, for example the excitation of the translational mode £j ' = | (j]^ 

sponds to the synchronous translation of the double kink whereas the excitation of the internal mode corresponds 
to the synchronous wobbling of both components. 



The meaning of the second mode is different. Indeed, the second equation (25b) is of the same type, so it can be 
solved by reducing it to the hypergeometric form (see, e.g. 
us introduce the ladder operators 



it _ 



d_ 
dx 



ntanh x; 



J). Here, however, we shall follow another route. Let 

(27) 



d_ 

dx 



n tanh i 



which allows us to transform the equation ( 25b ) to the form 



(a^a + n-2) & =u 2 & 



Here we introduced a new parameter 



1 1 /25 + 7k 

2 + 2 V 1 — k 



which is defined as a solution of the algebraic equation 

n(n + 1) 



6 + 2k 

1 - K 



Evidently, the ground state is defined as the state annihilated by the operator a, i.e. 

a£j = 0. 

Thus, up to a normalisation factor, the solution to this equation is 



AO) 

S2 



1 



cosh™ a 



(28) 



(29) 



(30) 



(31) 



(32) 



Here the integer n labels the internal modes of the double kink system, their number increases with coupling constant 
since 



K(n) = , 



n 2 + n — 6 
n 2 + n + 2 ' 



(33) 



Thus, for — 1 < k < there is only one bound state, 
k < 7/11 there are three internal modes, etc. Indeed 
U(x) = n(n + 1) tanh 2 



for < k < 3/7 there are two internal modes, for 3/7 < 
the properties of the modes are specified by the potential 
2 of the equation (|25bJ), it is getting deeper as the integer n increases. The number of 



bound states trapped by this potential then also increases. 



Note that the frequency of the lowest mode is w 



2, thus if n < 2, or equivalently, if k < 0, there is a negative 



mode in the spectrum and the lump kink becomes unstable with respect to perturbations. 

t (0) 



10 



we can see that the excitation of the lowest translational mode leads to decay of the double 
kink into two lump kinks which go off in contrary directions. On the other hand, the same perturbative distortion 



If the coupling constant 

k remains positive, the lump kink is stable and excitation of this mode = | — ^2°^) corresponds to the 

oscillations of the components of the double kink about the center of mass of the system. 

Our numerical studies confirmed this conclusion. An example of the decay of the double kink at K = —0.1 is 
displayed in Fig. 
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Figure 16: (Color online) Decay of the double kink at k — —0.1. Initial configuration is taken as <f>i(x,0) = utanh(x + 0.1) and 

4>2(x, 0) = wtanh(a; — 0.1). Excitation of the lowest translational mode £2 leads to decay of the double kink into two lump 
kinks moving in opposite directions. 



of the double kink at positive coupling k — 0.5 and consequent excitation of the wobbling collective mode does not 
destroy the configuration. 

The vibrational mode of the double kink has the following form (up to a normalisation factor): 



sinh x 
cosh™ x 



(34) 



The corresponding frequency depends on integer n as uj = y/Z{n — 1). Note, that action of the creation operator on 
the ground state (32) yields: 



2n sinh x 
cosh n+1 x 



(35) 



Evidently this is a (n + l)-th bound state. So in order to obtain the n-th bound state, the creation operator should 
act on the (n — l)-th mode. Indeed, the algebra of the ladder operators (27) is 



[a, a f l 



2n 



cosh 2 : 



(36) 



Therefore 



a) a + n — 2) a T = a T (a T a + n — 2 + 2n — In tanh^ x) = a) \ a) a + n — 2 + p 



;t - At (tf 



(37) 



where p = 2n + 1 and we introduce new ladder operators of degree (n — 1) (cf. Eq. (27)) 



(n — 1) tanha; , 



5 = - — h (n — 1) tanh x . 
ax 



(38) 
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Figure 17: Profiles of the lump kink components 4>i,cj>2 together with the profiles of the corresponding internal modes 771,772 
are presented at k = 0.5, the frequency of the internal mode is uj — 1.8537. 



(k) 

This, while the function ^ is an eigenfunction of the operator of degree n ( 28 1 and the corresponding eigenvalue 
is uj 2 , the action of the raising operator a transforms it to the eigenfunction of the operator of degree (n — 1) with 
eigenvalue lu 2 + p. 



B. Spectral structure of lump kink 



Investigation of the spectrum of the linear lump kinks fluctuations is a bit more complicated problem because, unlike 
in the double kink case, there is no analytical solutions for the classical solutions <j)\ and 4>2- Thus, our analysis relies 



on numerical methods, we have to find solutions to the system of equations (22) using shooting in two dimensions or 
applying the spectral Chebyshev method. 



Note the the system of the linearized equations (22 1 can be diagonalised on the spacial asymptotics where the fields 
approach the vacuum value: <j>i = (f> 2 = l/Vl — K as x — > ±oo. Introducing the same variables £1,^2 (24 1, we obtain 
then decoupled system of equations 



dx 2 '■ 4 )&=« 9 6. 



d , 1 + k 1 . 



(39a) 
(39b) 



Thus, there are 2 different excitations about the vacuum in our model (JlJ with 2 different masses mi — 2 and 

TO2 = 2y jz^) they are degenerated only in the special decoupled case k = 0. One of the interesting manifestation of 

the presence of two different excitations in the spectrum of linear fluctuations in our model is related with behavior 
of the lump kink under the influence of an incident wave. It is known in the usual one-component 4 theory the kink 
starts to accelerate in the direction of the incoming wave [33 with acceleration proportional to the 4-th power 
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Figure 18: Profiles of the lump kink components cj>i,(f>2 together with the profiles of the corresponding internal modes 771,772 
are presented at « = —0.5, the frequency of the internal mode is uj — 1.12245. 



of amplitude of incident wave. This effect is known as the negative radiation pressure. In our model \4>i\ — \<j>2\ 
symmetry restores exactly the 4 model so within this sector of theory the same phenomenon takes place. 

However in our system there is also a different possibility. Because of the difference in masses of small perturbation 
around the vacuum another mechanism of the negative radiation pressure (proportional to the square of the amplitude 
of the wave) can be observed [26] . In the case under consideration numerical simulations confirm that if the frequency 
of the incoming wave is u 2 = ra\ + k 2 the kink accelerates towards the source of the radiation, however if the frequency 
of the incoming wave is lo 2 = m 2 + k 2 the result is inverted, i.e. the positive radiation pressure is observed. 

Thus, the masses of the bound states of the lump kink are restricted by the lowest of the asymptotic masses mi and 
772,2- We found the solution of the eigenvalue problem numerically using the spectral Chebyshev method. The results 
are presented in Figs. |17|18| Evidently, the internal mode of the kink component is antisymmetric while the internal 
mode of the lump component is symmetric. Considering the dependency of the frequency of the internal mode on 
the coupling constant k we found that there is no second internal mode in the system, the second mass threshold is 
decreasing faster than the frequency of the oscillation mode (cf. Fig. 19 1. 



V. CONCLUSIONS 



Motivated by the recent interest in investigation of the remarkable resonance structures in the parameter space of 
the kink-antikink system, we have studied dynamical properties of the soliton solution of the two-component coupled 
</> 4 system. There are two different types of the solitons, the double kink and the lump kink, in both cases our numerical 
simulations of the collisions between the various types of the solitons reveal some chaotic resonance behavior which 
however, is rather different from the usual pattern of the KK collision in the one-component </> 4 model. 

Considering collision of two lump kinks in the same sector below some critical velocity we observed the resonance 
energy exchange between the sectors with a sequences of bouncing windows similar with the quasi-fractal dynamics 
observed in the usual <j> model. However the structure of the resonance windows in the former case is not so regular, 
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Figure 19: Frequency of the internal mode of the lump kink is presented as a function of coupling constant k. 

also some false windows appear. At the collisions of the lump kinks above the critical velocity we observed flip of 
the sectors and consequent scattering of the solitons. In the ultrarelativistic limit another process starts to dominate, 
the collision of the kinks in the first sector is accompanying by the collision of the lumps in the second sector which 
produces kink-antikink pairs. 

We also analysed the collision of the lump kinks in the different sectors and the most complicated process of the 
douple kink - lump kink collision. In the former case the resonance structure is very regular, it is qualitatively similar 
to the structure which appears in the one-component model. However the double kink - lump kink collision may 
lead to various results, there arc different channels of this process, the most interesting case is related with collisions 
at relatively high impact velocities, then a new type of self-similar fractal structure appears. These results need 
more rigorous investigation which is currently in progress, in particular it would be interesting to find an effective 
theory which could explain such a behavior. Evidently such a complicated behavior is related with the spectrum of 
excitations which may be excited in the collision of the solitons and affect the mechanism of the energy exchange 
between the solitons. 

Investigation of the linear stability of the soliton solutions showed that for certain range of values of the coupling 
constant, the double kinks are unstable with respect to linear perturbations, they decay into pair of two lump kinks. 

We also investigated the spectrum of perturbations of the soliton solutions on the two-component model. We found 
that in the case of the double kink configuration, in addition to the expectable counterparts of the translational and 
internal modes, there is a tower of bound states whose number depends on the strength of the coupling between 
the sectors. Evidently, excitation of these states in the process of the collision of the kinks will strongly affect the 
mechanism of the energy transfer. Considering the spectrum of linear fluctuation about the lump kinks we found 
there are excitations of 2 different masses, however there is only one internal mode in the system for all range of values 
of the coupling constant. 

It remains to systematically analyze the effect of interaction of the soliton solutions of the two-component model 
with an incoming wave in perturbation theory, it will answer the question if the effect of negative radiation pressure 
is also presented in this case. 

As a direction for future work, it would be interesting to construct an effective collective coordinate Lagrangian for 
the two-component model which will capture the most important degrees of freedom in the soliton collision. 

We thank Patrick Dorey and Wojtek Zakrzewski for enlightening discussions. Ya.S. is very grateful to Stephane 
Nonnenmacher for kind hospitality at the SPhT, CEA Saclay where part of this work was done. A.H. gratefully 
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Note that usual classical rescaling of the 2-component model in one spacial dimension does not allow us to absorb all the 
constants into rescaled scalar fields and couplings. 
[39] Note that it was shown recently [52] that the resonance windows in a model with a polynomial potential may appear due 
to existence of the collective bound states trapped by the kink-antikink pair. 



